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Abstract 

We present a first estimate of the cross-section for the exclusive process 
lHQDlKQl) ~^ P°lP°li which will be studied in the future high energy e + e~ — 
linear collider. As a first step, we calculate the Born order approximation of 
the amplitude for longitudinally polarized virtual photons and mesons, in the 
kinematical region s ^> —t, Q\,Q\- This process is completely calculable in the 
hard region Q\ , Q\ ^> Aq CD . We perform most of the steps in an analytical way. 
The resulting cross-section turns out to be large enough for this process to be 
measurable with foreseen luminosity and energy, for Q\ and Q\ in the range of 
a few GeV 2 . 
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1 Introduction 



The next generation of e + e~— colliders will offer a possibility of clean testing of QCD 
dynamics. By selecting events in which two vector mesons are produced with large ra- 
pidity gap, through scattering of two highly virtual photons, one is getting access to the 
kinematical regime in which the perturbative approach is justified. If additionally one 
selects the events with comparable photon virtualities, the perturbative Regge dynamics 
of QCD of the BFKL [1] type, should dominate with respect to the conventional partonic 
evolution of DGLAP [2] type. Apart from the study of the total cross section which has 
been proposed as a test of BFKL dynamics [3, 4], one can achieve similar goals by study- 
ing diffractive reactions. From this point of view the production of two J/^— mesons 
was studied in Ref. [5]. In this case the hard scale is supplied by the mass of the heavy 
quark. In the present paper we propose to study the electroproduction of two p— mesons 
in the 7*7* collisions. In this case the virtualities of the scattered photons play the role 
of the hard scales. As a first step in this direction we shall consider this process with 
longitudinally polarized photons and p— mesons, 

i*M 7lfe) - pdh) PL (k 2 ) , (1.1) 

for arbitrary values of t — {qi — ki) 2 , with s ^> —t. The choice of longitudinal polarizations 
of both the scattered photons and produced vector mesons is dictated by the fact that 
this configuration of the lowest twist-2 gives the dominant contribution in the powers of 
the hard scale Q 2 , when Q\ ~ Q\ ~ Q 2 . The measurable cross section is related to the 
amplitude of this process through the usual photon flux factors : 

da(e + e~ — > e + e~pp) 1 a 

ki^om = ofofSr 7/ 7 - pp) • (L2) 

where are the longitudinal momentum fractions of the bremsstrahlung photons with 
respect to the respective leptons and with P 7 / e (y) = 2(1 — y)/y for longitudinally po- 
larized photons. These double tagged events allow to access this hard cross-section [6]. 
According to our best knowledge this process has not been discussed up to now. At lower 
energy, some experimental data exist for Q\ small [7] and may be analysed [8] in terms 
of generalized distribution amplitudes [9] or transition distribution amplitudes [10]. 

In this paper we calculate the scattering amplitude of the process (1.1) in the Born 
approximation. In this way we get an estimate of the cross section and prove the feasibility 
of a dedicated experiment. Partial results have been presented in Ref. [11]. In the near 
future, we intend to extend this study by taking into account of BFKL evolution and 
transverse photon polarizations. 



2 Kinematics 

The process (1.1) at high energies can be visualized as in Fig.l. The physical picture of 
this impact representation is due to the presence of different time scales: a qq dipole is 
formed from the virtual photon, then interacts by exchanging t— channel gluons before 
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Figure 1: Amplitude for the process 7272 — > p\p\ m the i m P a ct representation. The 
blobs denote the vector meson distribution amplitudes. 



recombining in a meson. The Fig.l also explains the kinematics of the process (1.1). We 
parametrize the incoming photon momenta by introducing two light-like Sudakov vectors 
q[ and q' 2 related to the incoming particles, which satisfy 2q[ ■ q 2 — s = 2q± ■ q 2 . The usual 
s 7 » 7 * is related to the auxiliary useful variable s by s 7 . 7 » = s — Q\ — Q\. In this basis, the 
incoming photon momenta read 

/ Qi i 

Qi = Qi ?2 » 

s 

q 2 = q'2-%[- (2.1) 

s 



Their polarization vectors are 



Mi) + ^9l„> 



4 (2) = 7? + , (2-2) 

which is obtained after imposing the conditions = 1 and • e L ^ = 0. Because of e.m. 
gauge invariance, polarization vectors (2.2) can be effectively replaced by the second terms 
on the rhs of these equations. The momentum transfer in the t— channel is r = k\ — q±. 

We label the momentum of the quarks and antiquarks entering the meson wave func- 
tions as l\ and l[ for the upper part of the diagram and l 2 and 1' 2 for the lower part (see 
fig.l). 

In the basis (2.1), the vector meson momenta can be expanded in the form 

h = a(ki) q[ + — q' 2 + r± , 

a s 

2 

k 2 = (3(k 1 ) q' 2 + "j-^ - r ± . (2.3) 

Note that our convention is such that for any tranverse vector v± in Minkowski space, v 
denotes its euclidean form. In the following, we will treat the p meson as being massless. 
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a and (3 are very close to unity, and read 



«(*i) = 2l 1 --r 



1-4 



\ S (l-?)(l-f) 



i + 



V 



1-4 



\ s (i-¥)(i-f) 



(2.4) 



In this decomposition, it is straightforward to relate t = r 2 with r 2 = — r\. The 
corresponding relation is 



t 



Qi Q2 



2r 2 



1 



1 



1+ 1 



4r 2 



1 



I Q± ^ _ ^2 



- 1 



(2.5) 



or equivalently, 



r 2 = 



t + 



QiQl 



1 _ 5i) h _ 91 



1 



Q\Ql\ 



1 + 



—2— 



(2.6) 



From Eq.(2.5) the threshold for \t\ is given by \t\ min = Q\Q\js , corresponding to r± = 0. 
In the kinematical range we are interested in, the relation (2.6) can be approximated as 
r 2 = —t, in accordance with the usual result that r can be considered as purely transverse 
in the Regge limit. 

We write the Sudakov decomposition of the quarks entering the p mesons as 

h = z x q x + < ± i + z ± r± q 2 , 

Z\S 

/ x = z x q x - + Zir L = q 2 , 



, / . , ( l ±2 - z 2 r±) 2 , 
h = z 2 q 2 + I ±2 ~ z 2 r± q 1 



z 2 s 



- , j - H±2 - z 2 r±) 2 , 



z 2 s 



(2.7) 



where we have explicitly separated the transverse momenta l ±1 (I ±2) of quark and anti- 
quark forming the p^— mesons with respect to the momentum k\ (A^). In the following we 
shall apply the collinear approximation which consists in putting the relative momenta 
li± in Eq.(2.7) to zero at each qqp— meson vertex. The decomposition (2.7) is more easily 
understood in a slightly different Sudakov basis, which up to term linear in r± is obtained 
by the substitutions q[ — > q[ + r± and q' 2 — > q' 2 — r±, see Eq.(2.3) with a(k±) and (3(k 2 ) 
close to unity. In this new basis the p mesons have no transverse momenta and quarks 
transverse momenta have only relative momenta. 
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3 Impact representation 

The impact representation of the scattering amplitude for the process (1.1) has the form 
(see Fig.2) 




Figure 2: Amplitude for the process 71,72 — > p\p\ at Born order. The t-channel gluons 
are attached to the quark lines in all possible ways. 



M = is [ J^ q ^ p0 ^\k,r-k) J^ q2 ^ p °^(-k,-r + k), (3.1) 

J {2ix) A k (r — k) 2 

where J'^l(^)~*Pl ( kl ^(k, r—k) (J"fi(<)2)->PL ( k2 \k, r—k)) are the impact factors corresponding 
to the transition of 72(gi) — > p° L (ki) (72(^2) — ► Pl(^2)) via the t— channel exchange of 
two gluons. The impact factors are the s— channel discontinuities of the corresponding 
S— matrices describing the 7^ g — > p\ g processes projected on the longitudinal (nonsense) 
polarizations of the virtual gluons in the t— channel 

J^)^W{k,r_-k) = J ^s;l (qi) ^ L{kl) q 2 ^ , (3.2) 

jTiteHpiW^ -r + k) = J ^Stb^^qfq? . (3.3) 

Here the integration variables a and (5 are defined by the Sudakov decomposition of the 
gluonic momentum k 

k = aq[ + [3q' 2 + k± . (3.4) 

The amplitude (3.1) depends linearly on s, since these impact factors are s- independent. 
Calculations of the impact factors in the Born approximation are standard [12]. 3 They 
read 

J^ qi) ^ kl) (k,r-k) 
1 

f e 5 ab 

= J dzxzxzx^z^ifois—p^ -^-Qif p a(k 1 )P P (z 1 ,k,r,iJ 11 ) (3.5) 



3 Recently the forward impact factor of ~f* L {Q 2 ) — * p° L transition was calculated at the next-to- leading 
order accuracy in Rcf. [13]. 



5 



and 



1 

/e 5 ab 
dz 2 z 2 z 2 (p(z 2 )8n 2 a s ^—Q 2 f p p(k 2 )P P (z 2 ,k,r,p 2 ) , (3.6) 



where 



Pp(zi,k,r,Jil) ^ + ^ + -2 £2 + ^2 ( Zl r-k)2+^ (2 lI _k)2 +/i 2 ( 3 ' 7 ) 

is proportional to the impact factor of quark pair production from a longitudinally polar- 
ized photon, with two t-channel exchanged gluons. In formula (3.7) the collinear approx- 
imation l± — has been made. 

In this expression p\ = Q\ z\ z~i + m 2 and p 2 = Q\ z 2 z 2 + m 2 , where m is the quark 
mass. The limit m — > is regular and we will restrict ourselves to the light quark case, 
taking thus m — 0. Note that Eq.(3.6) can be obtained from Eq.(3.5) from the combination 
of two substitutions, firstly (zi, Q\) — ■> (-2 2 , Q2) and secondly (fc, r) — > (fc — r, — r). This 
second substitution is easy to understand since for the upper (lower) blob, the total 
t— channel outgoing momentum is r± (—r±), the outgoing gluon carries momentum k 
(k — r), and the incoming gluon carries momentum k — r (k). This substitution effectively 
corresponds to exchanging the third and the fourth term in Eq. (3.7), and leaves this 
impact factor invariant, due to its symmetry under z — > z, which is reminiscent of the 
Pomeron structure of the t-channel state (namely, the impact factor is even under C 
conjugation). Thus, only the first substitution is necessary to obtain the upper blob from 
the lower blob. 

In the formulae (3.5, 3.6), <fi is the distribution amplitude of the produced longitudi- 
nally polarized p— mesons. For the case with quark q of one flavour it is defined (see, e.g. 
[14]) by the matrix element of non-local, gauge invariant correlator of quark fields on the 
light-cone 

1 

(0\q(x) Y q(-x)\p L (p)) =f P lT J dze^-^^z) , (3.8) 



where the coupling constant is f p — 216 MeV and where the gauge links are omitted to 
simplify the notation. The amplitudes for production of p°'s are obtained by noting that 
|p°) = l/V2(\uu) - \dd)). 

The structures in the formulae (3.5, 3.6) arise due to the collinear approximation. 
Indeed, one can neglect the l± dependencies in the propagators, when computing the 
impact factors. Thus the only dependency with respect to this transverse momentum 
of the quark, in the Sudakov basis of the p— meson, is inside the wave function of these 
p— mesons. When integrating over the phase space of the quarks and antiquarks, this 
results in integrating the meson wave functions up to the factorization scale (which is of 
the order of the photon virtualities). We neglect here any evolution with respect to this 
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scale. Such integrated wave function is by definition the distribution amplitude of the 
p— mesons. For simplicity, we use the asymptotic distribution amplitude 



4>{z) = 6z(l-z). 



(3.9) 



The only remaining part the quark phase space integrations are the integration with 
respect to quark longitudinal fractions of meson momenta Z\ and z 2 . 
Combining Eqs.(3.1, 3.5, 3.6, 3.7), the amplitude can be expressed as 



M=is2n N ° ' 



i 



v . t a s a em a{ki) (3(k 2 ) f p Qi Q 2 / dz 1 dz 2 z x z x <p(zi) z 2 z 2 c/)(z 2 )M(z 1: z 2 ) , 
c o 



with 

M(z u z 2 ) 



(3.10) 



d 2 k 



k{r-k) 2 
1 



1 1 

+ 



x 



zjr 2 + /if zfr 2 + /if (zir — k) 2 + /if {z\r_ — k) 2 + /if 
1 1 1 



z 2 l 2 + a*! z \ r - + iA ( z 2L — k) 2 + iA ( z 2L — k) 2 + iA 



. (3.11) 



In term of this amplitude, the differential cross-section can be expressed in the large s 
limit (neglecting terms of order Q 2 / s) as 



It 



\M\ 2 

16 7TS 2 



(3.12) 



4 Cross section at Born order 



4.1 Forward case 

We begin with the simpler case t = t min (i.e. r = 0), where the final result for the function 
M(z\,z 2 ) (3.11) can be written in a rather simple form. The integral over k can readily 
be performed and gives 



M( Zl ,z 2 ) 



An 



Z\ Z\Q\ 

Z\ Z\Z 2 z 2 Q\ Q 2 (z\ z,\Q\ — z 2 z 2 Q 2 ) z 2 z 2 Q\ 



In 



(4-1) 



The amplitude M. given by Eq.(3.10) can then be computed analytically from Eq.(4.1) 
through double integration over z\ and z 2 . This is explained in appendix A.l. It results 
in 



M 



N 2 - 1 97T 2 

-is C N2 a 2 s a em a(k!) (5{k 2 ) f 2 

c 



6 R 



1 

R 



ln 2 R 



p 2 Q 2 Q 2 

+ 12 ^R-^lnR + 12 ( R + ^j + (s/? 2 + 2 (\n(l-R) In 2 R 

- In (R + 1) In 2 R - 2 Li 2 (-R) In R + 2 Li 2 (R) In R + 2 Li 3 (-R) - 2 Li 3 (R))] , 



(4.2) 
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where R = Q\/Q 2 . 

In the special case where Q = Q 1 = Q 2 , it simplifies immediately to 

N 2 -1 Q71- 2 
M tinin (Qi = Q 2 ) ~ -i 5 ^ a em tt(^i) m) ft (24 - 28 C(3)) . (4.3) 

The peculiar limits i? 3> 1 and i? 1 are of special physical interest, since they corre- 
spond to the kinematics typical for deep inelastic scattering on a photon target described 
through collinear approximation, i.e. the usual parton model. At moderate values of s, 
apart from diagrams with gluon exchange in t— channel considered here, one should also 
take into account of diagrams with quark exchange. We do not consider them here since 
we restrict ourselves to the asymptotical region of large s. 
In the limit R >> 1, the amplitude simplifies into 

tmin ~ a(W(* 2 ) /, - - j . (4.4) 

This result can be obtained directly by imposing from the very beginning the k± ordering 
typical of parton model. This is shown explicitly in Appendix A. 2. 



4.2 Non forward case 

In this section we will compute the amplitude (3.1), where we have been able to perform 
analytically the k±_ integrals. It involves the evaluation of a box diagram with two distinct 
massive propagators and two massless propagators (denoted I± ma m h below). We are not 
aware of any previous analytic calculation of such an integral. This gives us the possibility 
of studying various kinematical limits in variables Q\, Q 2 ,, t. 

M(zi, z 2 ) as defined in Eq.(3.11) is symmetric under z\ <-> Z\ and under z 2 <-> z 2 . 
Since on the other hand the distribution amplitude (f>(z) is also symmetric under z <-> z, 
M(zi, z 2 ) can be modified by adding any antisymmetric term, since the integration over 
Z\ and z 2 will automatically select its symmetric part. 

One thus writes 

1 

N 2 — 1 f 

M =is2n—j^ — a^aemaik^Pih) f 2 Q 1 Q 2 / dz x dz 2 z ± z x </>{zi) z 2 z 2 <p(z 2 ) M asym {zi, z 2 ) , 



(4.5) 



with 



M aS ym. (Zl, Z 2 ) — J 



d 2 k 



k(r-k) 



2 



1 1 



4r 2 + (il)(z 2 2 r 2 + n 2 ) (z 2 r 2 + ^ 2 )((z 2 r - k) 2 + /i 2 ) 

(4.6) 



1 1 

+ 



(z 2 2 r 2 + i4)((z lL - k) 2 + fx 2 ) {(z lL - k) 2 + fi 2 )((z 2 r- k) 2 + li 2 ) 
M asym .(zi, z 2 ) can be expressed in term of three kind of integrals, namely, 



8 



d d k 



'3m 



k{k-p) 2 

d d k 



k 2 {k-p) 2 {{k-a) 2 + m 2 ) 



and 



L 4m a m b 



d d k 



k it - p) 2 ((k - a) 2 + m 2 a )((k - 6)2 + m 2 



(4.7) 
(4.8) 

(4.9) 



where we use the dimensional regularization d = 2 + 2e. However, since we will effectively 
rely, for computing the various integrals involved, on a method which is applicable only 
for both UV and IR finite integrals, it is more efficient to rewrite M (3.10) in terms of 
another asymmetric two dimensional amplitude of the form 



M(zi, z 2 ) = - 



z 2 r 2 + p\ 



+ 



z\r 2 + n\ 



J^ 2 {z 2 ) - (1 <-> 2) 



where 



and 



k{k-L) 2 



1 



(fc — ra) 2 + /x 2 a 2 r 2 + /i^ 



+ (a <-> a) 



(4.10) 



(4.11) 



JAn lt i 2 {z\, Z2) — J 



d 2 k 



k(k-r) 



x 



((k - r_ Zl ) 2 + fi 2 Mk - r_z 2 ) 2 + /i 2 ) (z 2 r 2 + tf)(z 2 r 2 + fi 2 ) 



+ (z <-> z) 



(4.12) 



J 3fM and J4^ 1M2 are two dimensional integrals with respectively 3 propagators (1 massive) 
and 4 propagators (2 massive, with different masses). They are both IR and UV finite. 
Their computation by brute force technique, using Feynman parametrization, seems un- 
tractable in such a form (specially for J4 m//2 ). However, applying a trick inspired by 
conformal field theories, it is possible to compute these integrals. The basic idea is to 
perform special conformal inverse transformations, considered here in momemtum space. 
Although these two integrals, because of mass terms, are not conformal invariant, this is 
actually efficient after a suitable redefinition of the massive parameters. This is presented 
in the appendix A. 3. 

To complete the evaluation of the amplitude M., one needs to integrate over the quark 
momentum fractions in the p mesons z 1 and z 2 - In the general case, for arbitrary values 
of t, it seems not possible to perform the z\ and z 2 integrations analytically. We thus 
do them numerically. In course of them, we observe the absence of end-point singularity 
when 2:1(2) — > or 21(2) — > 1, since Pp as defined in Eq.(3.7) diverges like 1/z, 1/(1 — z) 
when z — > 0, 1. This leads to perfectly stable numerical integrations. 
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4.3 Results 



We use now the previous formulae in order to get prediction for production rate of diffrac- 
tive double p production. Running of as is a subleading effect with respect to our 
treatment. Anyway, we choose to replace a 2 s in the various formulae presented above 
by as(Qi) as(Q%) in order to fix the coupling, and use the three-loop running coupling 
a s (Q\) and a s (Q 2 2 ) with A^L = 305 MeV (see, e.g. [15]). 

In Fig. 3, we display the differential cross-section da^l'yl — > p° L p° L )/dt for vanishing 




2 4 6 8 10 

Ql/Ql 

Figure 3: Differential cross-section for the process 72 7£ — > p\p\ at Born order, at the 
threshold t = t m i n , as a function of Q\jQ\- The dots represents the value of the cross- 
section at the special point Qi = Q 2 , as given by the analytical formula (4.3). The 
asymptotical curves are valid for large Q^/Qf, as predicted by the asymptotical form 
(4.4) 

transverse t— channel momentum, i.e. t = t m i n , as a function of the ratio Q\/Q\. Curves 
are labelled by the values of Q\. The dots on the curves represent the values of the cross- 
section at the special point Q\ = Q2, which obviously correspond to the analytical formula 
(4.3). The cross-section dramatically decreases when Q\jQ\ increases at fixed Q\. For 
comparison, we show for each value of Q\ the asymptotical curve obtained by combining 
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Eq.(4.2) with Eq.(3.12). The complete result approaches quickly its asymptotical curve. 
However, in view of the strong decrease of the differential cross-section with increasing 
Q\lQ\i the asymptotical result seems to be of little interest for estimating data rates in 
the most favorable kinematics. 

The t-dependence of the differential cross-section da/dt is displayed in Fig. 4 for 
various values of Q = Q\ = Q2- 




2 4 6 8 10 

-t (GeV 2 ) 

Figure 4: Differential cross-section for the process 7^72 — > p\p\ at Born order, as a 
function of t, for different values of Q = Q\ = Qi. 



As anticipated, the cross-section is strongly peaked in the forward direction. This 
fact is less dangerous than for the real photon case since the virtual photon is not in the 
direction of the beam. However, the differential cross-section seems to be sufficient for 
the t— dependence to be measured up to a few GeV 2 . The comparison of the curves on 
Fig. 3 for t = t min with those on Fig. 4 leads to the conclusion that the phenomenological 
predictions obtained in the forward case will practically dictate the general trends of 
integrated cross-sections. 

Figure 5 shows the integrated over t cross-section as a function of Q 2 = Q\ = Q\. The 
magnitude of the cross-section seems to be sufficient for a detailed study to be performed 
at the linear collider presently under study. Note that we did not multiply by the virtual 
photon fluxes, which would amplify the dominance of smaller Q 2 . However, triggering 
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efficiency often increases substantially with Q 2 [4]. At this level of calculation there is 
no s-dependence of the cross-section. It will appear after taking into account of BFKL 
evolution. 



10 e 
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10 4 r 



10 3 t 



10' 



10 1 r 



10' 



2 . 5 



7.5 10 12.5 



15 17.5 20 

Q 2 = Qi = Q\ 



Figure 5: The integrated cross-section for the process Yl 7l 
function of Q\ = Q\. 



p L p L at Born order as a 



5 Conclusion 

This Born order study shows that the process 7* (<5i)t* {Ql) ~ > P°lPl can measured 
at foreseen e + — e _ — colliders for Q\ , Q\ up to a few GeV 2 . Indeed, a nominal integrated 
luminosity of lOOfb -1 should yield thousands of events per year, with Q 2 > 1 GeV 2 . This 
would open a new domain of investigation for diffractive processes in which practically all 
ingredients of the scattering amplitude are under control within a pertubative approach. 

In the near future we expect to include the transversely polarized photon contribu- 
tion, which should slightly enhance the non forward amplitude (this amplitude obviously 
vanishes at t — t min ). We also intend to incorporate the effect of BFKL evolution, since 
resummation effects are expected to give a net and visible enhancement of the cross sec- 
tion. 
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Finally, let us note that the elusive Odderon may also be looked for in 7*7* exclu- 
sive reactions [16], and that one may use the strategy developed in Ref.[17] to find it 
through its interference with Pomeron exchange which gives rise to charge asymmetries 
in 7*7*— > imiiTi. 



NOTE ADDED : After this paper has been accepted, two works [21] improve our analysis 
by taking into account of higher order effects. 
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A Appendices 

A.l Computation of the amplitude at t = t 



mm 



In this appendix we will prove the result (4.2) for the production amplitude at \t\ min = 
QlQl/s , still neglecting the p-mass. The amplitude (3.10) can be written as 

M = is \ 1 a 2 s a ern a(h) (3(k 2 ) f K , (A.l) 

where 

1 1 00 

K= I ' dz x ! dz 2 [ dk 2 — Z lT^ 2 t^t, (A.2) 

J J J (k 2 + z 1 z 1 Q 2 1 )(k 2 + z 2 z 2 Q 2 2 ) : 1 ; 





which reduces to 



z\ Z\z 2 z 2 ^ Z\ Z\Q\ 
Z\ Z\Q 1 — z 2 z 2 Q\ z 2 z 2 Q\ 



K = I dz x I dz 2 - ^ 2 _ 2 In - ^ . (A. 3) 



In the following, we denote R = Qi/Q 2 . We first reduce the computation of K to a 
one dimensional integral evaluation. Performing the change of variables X\ — 4 z\ Z\ and 
x 2 = 4z 2 z 2 , K reads 

K= ±- 2 f dXl dX2X2 (l \n(R 2X -±) . (A.4) 



After replacing the variable x\ by x = xi/x 2 , K reads 

K=jLi^= f X2 ~r== ( l ~ i — ~52~) M-R 2 ^) ■ (A.5) 
4 2 Q? J VT=T 2 J ^J\-xx 2 V 1 - R 2 xJ v ; K ' 
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One can now split the integration domain x 2 G [0, 1] xx G [0, l/x 2 ] as a; G [0, 1] xx 2 G [0, 1] 
and x G [1, oo] xx 2 e [0, 1/x] . The integral corresponding to the second domain is identical 
to the first one after exchanging Qi and Q 2 . Performing the integration over x 2 in the 
first domain, one gets for K 

l 

i r dx 

K = ^2 J ^72 h 6 v^(l + x) + (3 + 2x + 3x 2 )(2 ln(l + y/x) ln(l - x))] 
1 

1 ^ ■ 



x V 1 " 1^7^ J ln(jR } + (1 ~ 2) • (A ' 6) 

This one dimensional integral, after the change of variable t = y/x can finally be reex- 
pressed as 



K = 



i 

64^ / dt {^ { ' 2t + ln(1 + t] ~ ln(1 " t)} HRt) " 6 (1 + ^ 
o 



x ( - ^-i^) ln(i^) - | ln(l + f) In (i2f) + | ln(l - t) In 



Mi-«)M*«m . f fl A. (A7) 



1 + i?t /J V # 

The integrals arising from the two first lines of the previous expression, supplemented by 
the corresponding R — > 1/R contribution, are easily computed by integration by parts 
through logarithmic and polylogarithmic function Li 2 . Using Landen relation for Li 2 (see 
Chap.l of Ref.[18]), this simplifies into 



A= 1 



64 Q 1 Q 2 



6 [R+-\ki'R+ 12 I R--\]nR+ 12 [R+^ 



(A.8) 



The two last lines of Eq.(A.7) (denoted B in the following) contain terms of the generic 
form 

ln(a + bx) ln(c + ex) , , . , 

—gdx, (A.9) 



f + gx 

which can be reduced to the standard form (see Chap. 8 of Ref.[18]) 



/ 



In(l-y)ln(l-cy)^. (A. 10) 



This last integral is evaluated in Chap. 8 of Ref.[18], for a restricted domain in c. An 
analytic continuation of such a result for the whole complex domain in c can also be 
performed [19]. Let us first define 

<p(ct) = arg(e ia ) - a (A.ll) 
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for any real alpha. With this definition, 4^ is a winding number which counts the 
number of turns one has to make around in order to bring back a to a value inside the 
interval ] — Then, the results for the integral (A. 10) reads, for x real 

ln(l - y) ln(l ~ cy)j = Lis (j^f) + Li 3 (^) + ^(1) - Li 3 (l - cx) 

-y \ + ln(l - x) Li 2 (l - cx) - ln(l - cx) Li 2 (x) 

+ (ln(l - cx) - ln(l - x)) Li 2 + ^ ln(l - x) + \ lnc ln 2 (l - x) 

\ C / D ^ 

— « 7T ci In 2 c + % 2 7rci lnc ln(l — cx) — « 2 ttci lnc ln(l — x) — 4 7r 2 c 2 c 3 In c 

+ 2 2 7T (ci — c 4 — C5) ln(l — x) ln(l — cx) — 4 7r 2 c 2 C3 ln(l — x) + 4 7T 2 c 2 C3 ln(l — cx) 

— 7r (ci — c 4 — c 5 ) i ln 2 (l — cx) + « 7r (— ci + c 5 ) In 2 (1 — x) + i 4 7r 3 c 2 c 3 , (A. 12) 

where 

ci = 77- V? ( ar g( c - 1) - ar g( c ) - arg(l - x)) , c 2 = — p (- arg(c)) 

c 3 = 2^ ^ ( - ar s( 1 _ x )) ' Q = 2tt ^ ( ar s( c ) + arg(x)) 

and c 5 = tp (arg(x) + arg(c — 1) — arg(l — x)) . (A. 13) 

2n 

Ci, c 2 , C3 and c 4 take either or +1 values while C5 can take —1, or 1 values. From 
this formula one can get an analytical form for B. The obtained result is very lenghtly 
and contains a large bunch of In, Li 2 and Li3 of various rational functions of R. A first 
simplification occurs when using the Landen relation for Li 3 

/A f 2 1 i 

Li 3 (x) = Li 3 ( - J - — ln(-x) - - ln 3 (-x) + -ip (arg(l - x) - arg(-x)) ln 2 x , (A. 14) 

(see Chap. 6 of Ref. [18]), which effectively enables one to get combinations of Li 3 of the 
only arguments R and —R. 

Then, one reduces the terms containing Li 2 as combinations of Li 2 of the only argu- 
ments R and —R. This is possible after using Euler relations for Li 2 , which reads 

7T 2 

Li 2 (z) + Li 2 (l -z) = — -\nz ln(l - z) , (A.15) 



and 



TT 2 



Li 2 (,2) + Li 2 y-j = — --\n 2 z +i(arg(z- 1) - arg(l -z))\nz , (A. 16) 

as well as Landen relation 
Li 2 (^) + Li 2 (y^J = -^ln 2 (l-^) +i V (- arg(l - z)) (ln(l - z) - In z - i tt) . (A.17) 
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The most non trivial transformation is based on Hill formula (which enables one to expand 
the double variable function Li 2 (xy)) continuated in the whole complex plane (see [19] 
for details), 



Li 2 (xy) = Li 2 (x) + U 2 (y) - Li 2 



x(l - y) 



+ i<P (arg(l - y) - arg(l - xy)) jln 
+ i (p (arg(l — x) — arg(l — xy)) < In 



1 — xy 
1 — x 



1 — xy 
1-2/" 



- Li 2 
-In 
-In 



xy 



3/(1-3?) 
1 — xy 

1 — xy 

x(l - y) 
1 - xy 



-In 



1 — x 
1 — xy 



In 



i-y 

1 — xy 



(A.18) 



+ V (arg(l -y)- arg(l -xy))ip (arg(l - x) - arg(l - xy)) . 
After a long succession of painful simplifications, one finally gets 

B = - — i (3 J R 2 + 2 ) (ln(l -i?) In 2 i? - ln(i? + 1) In 2 i? - 2 Li 2 (-i?) lni? 

bAQiQ 2 \ R J 

+ 2U 2 {R) InR + 2Li 3 (- J R) - 2Li 3 (i?)) . (A.19) 

Using equations (A. 8) and (A.19) in order to express K = A + B and using Eq.(A.l), one 
proves the result (4.2). 

The limits R^> 1 and R <C 1 can be easily extracted using the asymptotical formulae 
for large argument of Li 2 and LI3. In the case of Li 2 they read 



7T 



1 



Li 2 (x) ~ — - In x + i(arg(x) - arg(-x)) lnx + — - ^— ^ - ^-^ + o ( — ) , (A.20) 



which reduces in the case of interest here to 

^ 1 1 



t • / \ 1 1 2 -i ^ 

Li 2 (x)~--ln x - in Inx + — — ■ 

2 6 x 4x z 

j • / \ l l2 • , 1 1 

Li 2 (x) ~ — -In x + 2 7r In x + 



9x 3 
1 



x 



Ax 2 9x 3 



x 



+ 



+ 



x- 



x- 



X" 



for x -> +00 , (A.21) 



for x -> -00 . (A.22) 



Similarly for Li 3 one gets from the functional relation (A. 14) the asymptotic expansion 



Li, (X) ~ -i - i ln(-x) + ± + ^ + ^ 

which reduces to 



+ 



x; 



T -/x 1, 1 .7T . 9 71-2 1 1 1 1 

L*(s)~--ln x-,-ln x+-mx + - + — + — 



fo I — 

X* 



and 
Li 3 (x) 



1, 3 2 A 1 1 1 
-In x+2-ln x+— InxH h— ^4 

6 2 3 a; 8x 2 



27x 3 



+0 I -3 



for x -> 00 , (A.23) 



for x — >• +00 , (A.24) 



for x -> -00 . (A.25) 



These asymptotical formulae immediately lead to 



M t 



+is- 



96tt 2 flnR 1 



R 6R 



(A.26) 
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A. 2 Computation of the amplitude at t = t m i n in the partonic 
approach 

In this appendix we rederive the asymptotic result (A. 26) in a way which corresponds to 
the usual parton collinear limit. In that limit, this simple results agrees with the fact that 
each loop of t-channel gluons gives rise to at most one logarithmic term. The logarithmic 
term corresponds to the leading logarithm approximation (LLA), while the constant term 
goes beyond this approximation. Let us show that this result can be easily obtained from 
the representation (3.7). Indeed, in the partonic approach, k\ is neglected with respect to 
any scale of the order of Q\, and symmetrically any scale of the order of Q\ is neglected 
with respect to k\. It is clear that in such an approximation, one immediately recovers 
the dominant contribution of (A. 26). The question arises how to define a precise 
prescription in order to get the full leading twist expression (A. 26). It turns out that the 
integral (A.2), provides the proper result 

1 f lnR 1 \ 
K ~ ~oFT~FT -5" - 7TB > ( A - 27 ) 



3QiQ 2 \ R GR 
after expanding the integrand K at leading order, namely 

c ~Fokv (A - 28) 

and then integrating the result from z 2 z 2 /R to Rz\Z\ . Let us justify this prescription. 
We rewrite the integral K of Eq.(A.2) as 

1 i 

2 2 OO 

K = -i- / dz, [ dz 2 [ du- z _lll^l (A.29) 

Q1Q2J J J (u + Rz 1 z 1 )(u + z 2 z 2 /R) : 1 ; 

000 

and separate the u = k 2 /(QiQ 2 ) integration as 



2929 

PRzizi 



du = du + du + du (A.30) 



0Rz 1 z 1 a^2£2 

where the parameters a and (3 are arbitrary. In the large R limit, f$Rz\Z\ < a^ 2 - for 
z \ < f ^j^r- Let us perform a systematic expansion of K in the limit where a and f3 satisfy 
R ^> a ^> 1 and R^> 1//3 ^> 1. We organize the expansion in such a form that the large 
R limit is taken first (which means the dominant twist approximation), and only then the 
large a and small (3 limit are taken. Decompose K = K 1 + K 2 + K 3 + K A + K 5 + K 6 , 
where 

1 a z 2*2 

2 9 ~W OO 

K\ = — — / dz 2 / dz\ \ du- —, 7— (A. 31) 

Q1Q2J J J (u + Rz lZl )(u + z 2 z 2 /R) 1 ; 

„ £2£2 
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1 a Z 2 Z 2 

2 ~W 



K 2 = _L_ / ^ 2 / / d«- Zi y Z2 _ (A.32) 

Q1Q2J J J (u + Rz 1 z 1 )(u + z 2 z 2 /R) 1 ; 



PRz-LZ-L 

1 a£2£2 

2 /3^^ PRz 1 Z 1 



K * = 7T7V I dz2 I dzi I du l — 1 p ZlZ \ Z2Z ^ — — — (A.33) 
Q1Q2 J J J (u + Rz 1 zi)(u + Z2Z2/R) 



00 

1. 1 
2 2 



K 4 = / dz 2 / dzi / du- — (A.34) 



a£2£2 0RZ1Z1 

R 2 



1 1 

2 2 



-^5 = „ „ ( dz 2 f dzi [ du- z ^2 z 2 (A.35) 



a Z 2 Z 2 



1 I n . Z 2 Z 2 

2 2 i? 



-^6 = ^ - / rf-22 / dz\ \ du- 1 Z } Z<1 2 ; — - . (A. 36) 

Q1Q2J J J (u + Rz 1 z 1 )(u + z 2 z 2 /R) 1 ; 

a Z 2 Z 2 

The integrals Ki, K 2 and K 3 are of order 1/R 3 and can thus be neglected at leading twist. 
It corresponds to the absence of end-point singularities in z variables and means that one 
could safely replace the lower bound of z\ integration in K4, K§ and K$ by 0. Let us 
focus now on the integral K 5 . The integration on u runs from a z 2 z 2 /R to (5 Rz 1 z±. Since 
Zi Z\ < 1/4 and z 2 z 2 < 1/4, in the limit » 1, a > 1 and /3< 1, one can safely expand 
the integrand of K 5 in powers of u/(Rz-\_z-\) and z 2 z 2 j{Ru). At leading twist, only the 
dominant term has to be kept, and K, can be approximated by (A. 28). Integration with 
respect to u, Z\ and z 2 then leads to 

K^(^.± + ± ln i). (A .37) 
5 \ R 6R 3R a J K J 

The logarithmic contribution In £ corresponds to boundary effect and is to be completely 
compensated by K 4 and K 6 , which behaves respectively as 1/Rln.p and 1/Rlna in the 
limit R >> 1, a >> 1 and C 1. This justify the assumption stated at the beginning of 
this paragraph. 



A. 3 Integrals 

In this appendix we collect all the generic integrals which appear in the computation of 
the Born amplitude. In addition to the integral I 2 defined in Eq.(4.7), we introduce the 
two following integrals with effective masses 



1 2m 



(A.38) 



kf((k-p) 2 + m 2 
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and 

d d k 



' hmamb J ((k-ay + mi)((k-by + mi)' (A ' 39) 

Let us first consider the case of the integral /2m • Using Feynman parametrization, one 
easily get 



1 2 X -l 



hm = tt 1+ T(1 - e)(p 2 + m 2 )- 1+t J daa^ 1 (l - c y^ 2 ) 



(A.40) 





In the massless case, one immediately obtains 

2tt 



/ 2 = ^(l + e(ln(V)-M/(l))). (A.41) 

For the non zero mass case, the a integration leads to the hypergeo metric function 2 Fi, 

^ = ^ f 1 - e, e, 1 + e; . (A.42) 

e{p z -\-m l y 6 \ p 2 + m z J 

After performing the Euler transformation z — > z/(z — 1) in the argument of the hyper- 
geometric function and then expanding the result in power of e, one gets 



'2m 



" ;i + e(ln7r-^(l)-lnm 2 + 21n(p 2 + m 2 ))). (A.43) 



e(p 2 + m 2 ) 



Let us now turn to the more general case where the propagators contain two (different) 
masses. In this case dimensional regularization is not necessary, and after straightforward 
calculations, using Feynman parametrization, one obtains 

J ^ (A 44) 

- h) 2 + (m a - m b ) 2 )((a- b) 2 + (m a + m b ) 2 ) 



x In 



(a - b) 2 + m 2 a + ml + yjlja - b) 2 + (m a - m b ) 2 )((a - b) 2 + (m a + m b ) 



(a - b) 2 + ml + ml— a/ ((a - b) 2 + (m a - m b ) 2 )((a - b) 2 + (m a + m b ) 2 ) 



For the purpose of our computation, we will need the previous integral only for the special 
case where a and b are collinear. 

In that case, one obtains, with r = |r|, 



J: 



2af3 



d d k 



((k - ra) 2 + a 2 )((k - rb) 2 + (3 2 ) 



7T r\a - b) 2 + a 2 + (3 2 + V\ 
~ V\ 11 r 2 [a - b) 2 + a 2 + (3 2 - V\ ' 

where we introduce the notation 

A(x, y, z) = x 2 + y 2 + ^ 2 - 2xy - 2xz - 2yz , (A.45) 
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which enables us to define, for the purpose of our computation, 



A = A(-r 2 (a - b)\ a 2 , (3 2 ) = {a 2 - (3 2 ) 2 + 2{a 2 + (3 2 )r 2 {a - b) 2 + r 4 (a - bf . (A.46) 
Let us first consider the integral I 3m . We start from the following identity 

r d 2 k ( 1 1 1 l \ 

J k^(k—p) 2 \(k — a) 2 + m 2 a 2 + m 2 (p — a) 2 + m 2 (k — p + a) 2 + m 2 J 



_i i \ r d d k 

a 2 + m 2 (p- a) 2 + m 2 ) J k 2 (k- p) 2 



d d k f d d k 

+ / To-. = TT + 



J k 2 (k-p) 2 ((k-p + a) 2 + m 2 )' [ ' ' 



kf(k-p) 2 ((k-a) 2 + m 2 ) J k 2 (k-p) 2 ((k-p + a) 2 + m 2 ) 

This identity relates a finite expression on the lhs with a sum of dimensionally regularized 
integrals. After shifting integration variable in the last integral on rhs of Eq.(A.47), one 
obtains 



/ 



d d k 1/1 1 \ f d d k , kin . 

+ 7Z ZV2 , ^ \ TTTZ^T, (A.48) 



k 2 (k-p) 2 ((k-a) 2 + m 2 ) 2\a 2 + m 2 (p-a) 2 + m 2 JJk 2 (k-p) 2 



i r d 2 k / l j. i i 



2 J k?(k — p) 2 \(k — a) 2 + m 2 a 2 + m 2 (p — a) 2 + m 2 (k—p + a) 2 + m 2 

which expresses I^m in terms of the divergent integral I2 already calculated and of a 
finite integral whose computation is our next task. The method which we use is the 
generalization for the massive case, in momentum space, of the technique of calculation 
of massless two dimensional diagrams in the coordinate space encountered in conformal 
field theories [20]. 

The essential point is to perform in two dimensional finite integrals a conformal trans- 
formation of type l—*l/l on the integration variables and vector parameters, and of type 
m 2 — >• 1/ m 2 for the dimensionful parameters. This transformation reduces the number of 
propagators. 

Let us illustrate this method in the special case where a and p are collinear, which is 
of practical interest for our computation. We will thus focus on the finite integral in rhs 
of Eq.(A.48), which turns out to be J 3m , as defined in Eq.(4.11). The transformation 



K R 1 

k -> r -> m—>— (A.49) 

- K ~ R 2 ' M 



reduces the number of propagators and gives 

f d 2 k f 1 1__ 

3m ~ J k 2 (k - r) 2 \(k- ra) 2 + m 2 rHm 2 ^"" ° j 

n2 f d 2 K ( K 2 R 2 1 

= R / ti? ~ D N2 j^w - -^r, — 2 + ° ^ ° • A - 50 

7 {K-R) 2 \(R-aK) 2 + ^fi- a 2 r 2 + m 2 
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After performing the shift of variable K_ = R + k/ and then finally making the inverse 
transformation 

fc'^A r^Z M-,-, (A.51) 

k rf m 

we end up with 



1 fd 2 k 



hm ~ r 2 J k 2 



{r + kf 1 
_ " +(a^a) 



a 2 r 2 + m 2 



I m^\(( h r r 2 aa-m 2 \ , m 2 r 4 

(A.52) 

The computation of this integral can now be performed using standard Feynman param- 
eter technique. It results in 

2tt / / 1 1 \ , r 2 a 2 + m 2 

•JZrn — — T } I o 2 i 2 2-2 i 2 J 2~2 i 2 (A.OOJ 

(I 1 2 \ (r 2 a 2 + m 2 )(r 2 a 2 + m 2 )\ , A . 

+ 2 2^ 2 + 2-2^ 2 + — 2 ln " ^4 " \ " A " 54 

Let us now focus on the /4 mm integral. We start with an analogous identity as those 
in Eq.(A.47), which leads to 



d d k 



t it - P) 2 ((k - a) 2 + m 2 a )(k - b) 2 + m\ 



(A.55) 



i 1 \ [ d d k 1 

771 _.\9 i 779x777. t\9 i -_9.\ ) / t 2/7 To ^ ' 



2 V(a 2 + m2)(6 2 + m2) ((p - a) 2 + m 2 a )((p - b) 2 + m 2 ) ) J k 2 {k-p) 2 2 
where 

f d 2 k t 1 1 

"7 k 2 (k-p) 2 \{{k-a) 2 + m 2 a ){{k-b) 2 + m 2 h ) (a 2 + m 2 ) (b 2 + m 2 ) 

1 1 

+ 



((p-a) 2 + m 2 )((p-6) 2 + m 2 ) ((fc - p + a) 2 + m 2 )((A; - p + 6) 2 + m 2 ) 

This identity express the integral l^ ma m b as a sum of a IR divergent integral which has 
already been computed in Eq.(A.41) and of a finite integral J, on which we will apply the 
same trick based on conformal transformations. Once more, although our method can 
be applied to the previous integral J, we will restrict ourselves to the more simple case 
where a and b are collinear. It thus means that we need to compute the integral as 
defined in Eq.(4.12). We apply on three successive transformations accompanied 

by appropriate changes of variables as described above: conformal transformation, shift 
of integration variable and inverse conformal transformation. This gives 

j 1 f d d k 
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\\dL — U ) 

<r»fc + rf) ((t + £^I) 2 + p^,) (r»% + A) ((i + r^gf) * + 

"(^ 1 + / ,;)( r ^ 2 + ^) + (z " S) J • 

In this expression, the integral which remains to be computed has the form 

T - f d 2 k((k-r_) 2 ) 2 

J ^ J k 2 [(k - ra) 2 + a 2 ] [(k - rb) 2 + (3 2 } ' 1 j 

This integral can be rewritten in the form 

J*# = JZ + J Z (A.59) 

where 

7^ = f d kk 

3a/3 7 [(£ - La? + a 2 ] [(fc - rb) 2 + f3 2 } 1 ' ' 

and 

7 « / ^[((fc-r) 2 ) 2 -(A; 2 ) 2 ] 

30/3 y A; 2 [(A;-ra) 2 + a 2 ] [(A; - r&) 2 + /3 2 ] ' l ' ' 

Jfajs is IR finite but UV divergent. On the contrary Jj^ is UV finite but IR divergent. 
The calculation of J^p is standard, although somehow lenghty, and leads to the result 

V\ ^ r 2 (a-b) 2 + a 2 + (3 2 + V\ 
2r 2 (a-b) 2 n r 2 (a - b) 2 + a 2 + (3 2 - V\ 
[V\-a 2 + P 2 - r\a 2 - b 2 )] 2 V\ - r 2 (a - bf - a 2 + (3 2 
4r 2 {a-b) 2 V\ ' n y/X + r 2 {a - b) 2 - a 2 + f3 2 

[V\-P 2 + a 2 - r\b 2 - a 2 )] 2 ^ y^A - r 2 (q - b) 2 + a 2 - p 2 \ 
Ar 2 (a-b) 2 V\ 11 V\ + r 2 (a-b) 2 + a 2 - (3 2 j ' 

The evaluation of can be simplified by expressing it in the following way: 

jir = 2 A +r 4 f tk 

3«/3 J -^k-ra) 2 + a 2 ][(k-rb) 2 + f3 2 } J k 2 [(k-ra) 2 + a 2 }[(k-rb) 2 + (3 2 } 

2 C d d kk ■ r T d d kk ■ r 

" r J k 2 [(k -ra) 2 + a 2 } \{k - rb) 2 + (3 2 } ~ J [{k - r_a) 2 + a 2 } [{k - rb) 2 + (3 2 } 

+ A f ^Cfe-r) a (A64) 

J k 2 [(k-ra) 2 + a 2 ][(k-rb) 2 + P 2 ]' V ' ; 
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Each of these integrals is regular, except the second one, which diverges like 



jIRdiv. _ 
J 3o/3 — 



e r(l + e) (r 2 a 2 + a 2 ) (r 2 b 2 + (3 2 ) ' (A ' 65) 

The evaluation of each of these integrals is lenghty but straighforward after using Feynman 
parametrization. 

Since IR and UV divergences cancel when evaluating the integral I defined in Eq.(A.57), 
we effectively need the finite part of J3 a /3, defined as the remaining term after removing 
the pole in 1/e, which sums up the finite part of Jj^ and of J^. With this definition, 
one can express J^ xm as 

J 4miM2 = (A.66) 



2 (r 2 z 2 + fij)(r 2 z 2 + /i 2 .) Jia(B { z\r 2 + /i 2 " ' z 2 r 2 + /i 2 " ' r 2 z\ + /i 2 ' r 2 z 2 + /x 2 ' 



finite 



f -z 1 z 1 r 2 + ji\ -z 2 z 2 r 2 + n\ 



r 2 \x 2 



jfinite 



( —z 1 zir 2 + n\ -z 2 z 2 r 2 + fx 2 . 



r 2 /ii 



r 2 [x 2 



+ r 2 (r 2 z 2 + fi 2 )(r 2 z 2 + fi 2 ) 3a/3 V z 2 r 2 + fi 2 ' z 2 r 2 + fi 2 ' r 2 z 2 + /i 2 ' r 2 z 2 + /i 2 ' 
with 



J 3a™ <e ( a > & > & r ) 



= 7T 



(a-p) 2 (a + py 
(a — 6) 2 r 2 
1 



+ 



In 



(A.67) 

4^—|^ - 2(a 2 + /? 2 ) - ((a - 6) 2 + 4(a + b) - 12)r 2 
a — o 

r 2 (a-b) 2 + a 2 + (3 2 + V\ 



x 



r 2 a 2 + a 2 r 2 b 2 + (3 2 
r 2 (a-b) 2 + 2(a 2 + (3 2 )-2 
r 2 b 



r 2 (a - b) 2 + a 2 + (3 2 - V\ 
a(3 2 -ba 2 



+ 2abr 2 + r 2 



a 2 -[3 2 + (a 2 - b 2 )r 2 



r 2 a 



+ 



r 2 a 2 + a 2 r 2 b 2 + (3 2 



-4 +2- 



(a 



a-b ' ' ' b(r 2 a 2 + a 2 ) - a(r 2 b 2 + (3 2 ) 

2 !3 2 ) 2 + 2{a 2 + l3 2 ) 2 r 2 {a-b) 2 + r\a-bY 1 



b (r 2 a 2 + a 2 )-a (r 2 b 2 + (3 2 ) 



a — b 



(a 2 -/? 2 ) 2 
r 2 (a — b) 2 



In 



v / A + r 2 (a-6) 2 + a 2 -/? 2 



+ 



a 2 /? 2 - 6 2 a 2 + a&^A 



r 2 a 



V\-r 2 (a-b) 2 + a 2 - (3 2 

r 2 b \V\ + a 2 -p 2 -r 2 (b 2 -a 2 ) 2 



ab(a - b)'a(r 2 b 2 + (3 2 ) - b{r 2 a 2 + a 2 ) 

2 -b(3 2 



2y/X \r 2 a 2 + a 2 r 2 b 2 + (3 2 J a(r 2 b 2 + (3 2 ) - b(r 2 a 2 + a 2 ) a-b 



{a 2 - l3 2 ) 2 (a 2 + b 2 )r 2 



2r 2 (a-b) 2 V\ 
x In — -— + 



2v/A 
1 



a + 6 
2(a - 6) 



a(r 2 !> 2 +,S 2 )-!>(r 2 a 2 + a 2 ) 



2 

r a 



ab 2(r 2 a 2 + a 2 )(r 2 b 2 + (3 2 ) 2 
r 2 b 



+ 



r 2 a? + a 2 r 2 b 2 + (3 2 



-4 + 



a 2 r 2 + a 2 



b 2 r 2 + (3 1 



In 



r 2 a 2 + a 2 



[r 2 a 2 + a 2 )(r 2 6 2 + (3 2 ) - abr 4 (r 2 a 2 + a 2 )(r 2 b 2 + (3 2 ) 



r 2 b 2 + (3 2 ab(r 2 a 2 + a 2 )(r 2 b 2 + f3 2 ) 



where appropriate additional lnr 2 terms have been introduced after extracting the fi- 
nite part of J^i and Jj^ in order to write the final result as made of logarithms of 
dimensionless quantities. This is possible since the final result Ja Ui u 2 is UV and IR finite. 
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